R parity violating SUSY explanation for the CDF Wjj excess 
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Recently the CDF Collaboration has reported a statistically significant excess in the distribution 
of the dijet invariant mass between 120 — 160 GeV in Wjj event sample in 4.3 fb -1 of data and 
later confirmed with 7.3fb _1 of data, which has generated considerable interest. We offer a possible 
explanation of this observation in the general framework of MSSM with imparity violation through 
resonance production of v r decaying into the LSP f\ and W boson. We also give the predictions of 
this scenario for the LHC operating at 7 TeV center of mass energy. 
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The CDF Collaboration at the Fermilab has recently 
reported a 4.1ct excess in the dijet invariant mass dis- 
tribution (120 GeV < Mjj < 160 GeV) in exclusive 2 
jets + lcpton + missing transverse energy ($t) events 
in 7.3fb _1 of data from pp collisions at 1.96 TeV where 
the lepton (I = e, u) and Jfix arc compatible with the 
decay of a real W PJ. A Gaussian fit to the excess re- 
gion of Mjj indicates the peak to be around 147 GeV 
and it is wider than the Z boson of the Standard Model 
(SM). There are few hundred events in the excess region 
of Mjj, which correspond to cr{Wjj) ~ O (pb). Such 
a large cross-section immediately rules out the possibil- 
ity of the indication of the SM Higgs boson H, because 
a(pp ->• W ± H) x BR(H -> 66) ~ 12 fb. The CDF anal- 
ysis also reported no significant deviation from the SM 
expectation of Mjj distribution in Z+jets and bb£~ $t 
samples. This result is in agreement with the earlier pub- 
lished one [3 by the CDF Collaboration with 4.3 ftT 1 of 
data sample, increasing the significance of the previous 
result. The position of the Mjj peak and the expected 
cross section for the process has not changed much. 

This anomalous excess in the Wjj event sample has 
naturally generated a lot of interest in the particle physics 
community, as this could be the long-anticipated di- 
rect signal of new physics beyond the SM. However, the 
D0 Collaboration, using their 4.3 fb -1 of data sample, 
claimed that the distribution of Mjj in such events is 
consistent with the SM prediction Q. Nevertheless, one 
should be rather cautious in either believing or disbe- 
lieving the claims made by the CDF or the D0 Collabo- 
rations. Before coming to any concrete conclusion, one 
must take into account the different methodologies used 
by the CDF and the D0 Collaborations to analyze their 
data samples, estimations of different systematics in the 
SM backgrounds, especially the QCD background. The 
D0 Collaboration by simulating pp — ► WH — > tvbb pro- 
cess to model acceptance and efficiency, put an upper 
limit of 1.9 pb on the cross-section of anomalous dijet 
production at 95% CL for Mjj = 145 GeV Q. As a result 
of this, their analysis does not rule out the possibility of 
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new physics interpretation of the CDF dijet excess with a 
cross-section less than 1.9 pb [4|]. On the other hand, one 
can argue that the CDF excess in dijet events could be 
due to some incorrect modeling of the QCD backgrounds 
in that mass window. To resolve this issue a joint task 
force has been formed [5j and till we reach a definitive 
settlement of the issue it might be interesting to explore 
the different theoretical avenues that could account for 
the dijet anomaly. Many suchpossibilities have been al- 
ready discussed recently 0, [a-tlO] . 

To explain the intriguing CDF result, we propose 
a scenario in i?-parity violating (RPV) supersymmet- 
ric (SUSY) model, where the lightest stau (fi) is the 
lightest supersymmetric particle (LSP) and the next-to- 
lightest SUSY particle (NLSP) is the corresponding sneu- 
trino, v T . All the other superpartners are rather heavy 
(O (TeV)). The v T may be produced resonantly via im- 
parity violating coupling and subsequently decays into a 
real f\ and an off-shell W through i?-parity conserving 
charged-current interaction. The real f then decays into 
pair of jets via the same i?-parity violating coupling and 
the off-shell W decays into tv. Such resonant sleptonpro- 
duction and decay has been briefly discussed in Ref. [7fl in 
the context of CDF Wjj anomaly with a smaller mass 
splitting between the sneutrino and the stau (only due to 
the D-term splitting). As a result of this, the lepton from 
the decay of the virtual W is rather soft to satisfy the 
CDF criteria of lepton selection which arises from the de- 
cay of real W. In another analysis [9(, authors assumed 
that the lepton and the neutrino come from the decay 
of the LSP charged slepton involving i?-parity violating 
couplings and not from the W-boson. The slcptons are 
produced in the decay of the pair produced neutral winos. 
In addition they assumed that the wino and the charged 
slepton are nearly degenerate in mass and showed that 
this scenario can explain the CDF dijet anomaly. 

In our analysis of the proposed signal pp — >• v T — >• 
W~f^~ — > t~vijj, we consider a large mass splitting be- 
tween v T and f i . This may be achieved if one considers 
the fi to be a mixture of left-chiral stau (tl) and right- 
chiral stau (tr) and hence the mass splitting between the 
v r and T\ can be made significantly larger (~ 200 GeV or 
so). To have (Mjj) in the correct ball-park, we consider 
Mfj in the range 140— 150 GeV. Now, the sneutrino {p T ) 
produced through the i?-parity violating coupling X' 3 n 



decays into fi and an on-shell W, which then eventu- 
ally lead to a pair of hard jets (from fi decay) and a 
charged lepton and neutrino (from real W decay) in the 
final state. 

In our numerical analysis, we use the CTEQ6L parton 
distribution function [ll[ with factorization scale Q = 
y/I/2, where VI is the parton level CM energy. We also 
smear lepton and jet energies according to 



cr(Ej) _ at 

Ei ~ y/E/GeV 



75% and b,, = 3% \l 
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Ej ^E/GeV 

where, at = 13.5%, be = 2%, 

After smearing we apply the selection criteria used by 

the CDF Collaboration, viz. 

• lepton (i = e, fi) and Et due to the neutrino 
should be consistent with decay of a W boson: 
E|(E£) > 20 GeV, \rj e \ < 1.0 and E T > 25 GeV 
with transverse mass Mt(I^) > 30 GeV. 

• Two jets with E J T > 30 GeV, cone size AR = 0.4; 



|?y J | < 2.4 such that 
40 GeV. 



Vh ~Vh l< 2 - 5 and (pt)jj > 



• | A0(Er,ji) |> 0.4 so that Et should not be due 
to mismeasurcment of jet energy 

• lepton and the jets should be isolated, AR(j,^) > 
0.52. 

With all these, we are now well equipped to explain 
the observations of CDF. We set M Pt = 300 GeV, M fl = 
145 GeV. From this mass spectrum, one can see that 
v T can be produced in s channel resonance via the X' 3ll 
coupling 1 , which then decay into the LSP fi and the 
gauge boson W. The LSP fi then decays into a pair of 
jets via the same R parity violating coupling and the W 
boson decays semileptonically, which eventually leads to 
jjivt final state with low invariant mass ( Mjj sa 140—150 
GeV). 

It has been shown in [jj that the bound on X' 311 cou- 
pling from dijet resonance searches both at the CDF and 
UA2 experiments allows one to take a value of this cou- 
pling to be < 0.3. In our analysis we take the value 
to be 0.1. With these set of parameters, we estimate 
&~lo(pp — ^ v T -> W~t^ — > jjW) = 1.296 pb, which is 
consistent with the 95% CL upper limit on a(Wjj) by the 
D0 Collaboration [3|]. This signal cross-section (1.296 pb) 
eventually comes down to 71 fb after taking into account 
the leptonic (e & fi) branching ratio of W boson and im- 
posing all the selection cuts mentioned above. There is 
no branching ratio suppression from f\ — > jj decay, as f\ 
decays into a pair of jets with 100% probability, which 



1 Here, we assume only one i?-parity violating coupling is domi- 
nant at a time. 
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FIG. 1: The top plot shows the distribution of dijet invari- 
ant mass (Mjj) after the CDF selection criteria have been 
applied: the points with error bars are taken from CDF data 
(Cint — 7.3 fb' 1 ), the histogram (red) shows our proposed 
signal events and a Gaussian is fitted to the histogram (blue 
line). The bottom plot shows distribution of Mwjj from our 
signal events. 



we think is a rather good approximation given the fact 
that all other SUSY particles are heavy. 

This cross-section corresponds to 518 (305) events with 
dnt = 7.3 fb _1 (4.3 fb -1 ). The M M distribution accord- 
ing to our proposed scenario is shown in Fig. [TJ (top 
plot). It is very clear that the signal histogram agrees 
reasonably well with the CDF data points. It may be 



noted that detailed implementation of detector effects is 
beyond the scope of this letter. Hence, instead of quan- 
titative comparison of the number obtained in this anal- 
ysis, which is somewhat larger than that of the CDF, 
we restrict ourselves to a qualitative comparison of the 
features of the distributions obtained. In the top plot 
of FigfTJ we normalize the signal histogram by the ratio 
of CDF peak value over the signal peak value. We also 
fit the signal histogram by a Gaussian and the mean of 
the Gaussian is given by My = 145 GeV. The CDF Col- 
laboration has given plots for several kinematic variables 
from their 7.3 ftT 1 data 0. In the bottom plot of Fig HI 
we display the invariant mass distribution of Mf„jj sys- 
tem. The longitudinal momentum (p z ) of the neutrino is 
obtained by constraining (pi + p„) 2 to M^ and extract- 
ing the two possible solutions. To construct the M^„jj 
we take the smallest solution between the two. The CDF 
Collaboration has shown a plot of Mf„jj after subtracting 
the SM background [13|. The plot shows excess in the 
region of our interest (250 GeV - 300 GeV). However, 
in their earlier analysis with 4.3 fb -1 data sample, this 
Mf„jj distribution was compatible in shape with the SM 
background only hypothesis Q. As a result of this, one 
requires more statistics to reach a definite conclusion in 
M^yjj distribution. 

At this stage it is also important to look at some other 
interesting predictions of this scenario in the context of 
Tevatron. For example, the pair production of f\ and 
subsequent decay of individual fi via X' 311 coupling can 
lead to four-jet final state. We find that the leading or- 
der cross-section for 145 GeV fi pair production is of 
the order of (~ 4 — 5 fb), which is rather small to be 
seen at the Tevatron with the available statistics. The 
fi can also be produced in association with Z and Higgs 
bosons H ± ,h°,H°,A° via Ajj X1 coupling. As before, fi 
decays into a pair of jets via the same i?-parity violat- 
ing coupling, while Z and Higgs bosons can have sev- 
eral possible decay modes, which eventually lead to some 
spectacular signatures. Once again, for our scenario, 
the production cross-sections for these processes are too 
small to have any relevance at the Tevatron energy. For 
example, we find that the leading order cross-section 
<jlo{pP — > t\ + Z) ~ 10~ 2 pb and we expect that the 
cross-section for a(pp -4 fi$°), (<f> = H ± , h°, H°, A ) 
would be of the same order or even smaller (for heav- 
ier Higgs bosons) than 10~ 2 pb. In addition to these, 
one can have associated production of a v T and fi and 
this can lead to the W + 4j signal where two pairs of 
jets will show peaks in their invariant mass distributions. 
Note that in this case the 2 — > 2 production process is 
not suppressed by the small i?-parity violating coupling. 
Similarly, pair-produced v T $> can lead to WW + 4j sig- 
nal at the Tevatron. Note that, with the mass spectrum 
considered here to explain the CDF dijet anomaly, one 
should not expect any significant number of events from 
those two processes with the present luminosity at the 
Tevatron. 

In conclusion in this letter we have analyzed the re- 



cently reported dijet invariant mass excess at 4.1er in 
M j:j ~ 120-160 GeV by the CDF Collaboration in 
Wjj events with 7.3 fb - data. We have proposed an 
i?-parity violating scenario with X' 3ll as the dominant 
coupling, where v T and fi are the NLSP and the LSP 
respectively. In this scenario, the mass splitting between 
the NLSP and the LSP is of the order of hundred GeV, 
such that NLSP can decay into the LSP and a real W 
boson. The resonant production of the NLSP can lead 
to a final state ivjj where, the pair of jets coming from 
the i?-parity violating decay of the LSP, show a peak at 
Mjj ~ 145 GeV with a rather good agreement with the 
dijet invariant mass distribution as shown by the CDF 
Collaboration at 7.3 fb - data. We have also found an 
s channel resonance at Mg V jj ~ 300 GeV, corresponding 
to our chosen value of the NLSP mass. This result also 
appears to be very close to the CDF plot of Mg v ^ after 
subtracting the SM background. However, for a defini- 
tive conclusion one requires higher luminosity. We note 
in passing that apart from the possible excess observed in 
Mt vj j distribution around 250 GeV - 300 GeV, the CDF 
excess in Wjj events in 4.3 fb" 1 as well as 7.3 fb -1 data 
samples are almost consistent with each other. We have 
also discussed other possible signatures of this scenario 
at the Tevatron. 

Let us also comment briefly on the implications of this 
scenario at the LHC, which is accumulating pp collision 
data at y/s = 7 TeV. We have found that the produc- 
tion cross-section for Wjj from pp — > v T — >• t\ + W is 7.3 
pb. At the LHC, the major SM background processes for 
this signal comes from the W + n jets process (n > 2). 
In the wake of the CDF result the ATLAS Collaboration 
performed a similar analysis with data from pp collisions 
at Vs = 7 TeV Q. They have analyzed 33 pb -1 of data 
collected in 2010 and did not observe any significant dif- 
ference between SM expectation and observation in the 
mass range of interest. The estimated W + n jets back- 
ground is approximately 20 times higher than the rate 
measured at the Tevatron [14| . After applying the same 
set of cuts as discussed before, the signal and the SM 
background cross-sections are 1.2 pb and 34 pb respec- 
tively. As a result of this, the signal significance of our 
scenario is < la at 33 pb~ of data. We know that both 
CMS and ATLAS have already collected L mt > 1 fb -1 
of data per experiment and are likely to collect a few 
ftV 1 of data by 2012. Hence we may expect that with 
the increase of luminosity it might be possible to have 
a definitive conclusion on this issue from both CMS and 
ATLAS Collaborations. 
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